AN INVERSE PROBLEM FROM SUB-RIEMANNIAN GEOMETRY 



THOMAS A. IVEY 
To the memory of Robert B. Gardner 

Abstract. The geodesies for a sub-Riemannian metric on a three-dimensional contact 
manifold M form a 1-parameter family of curves along each contact direction. However, a 
collection of such contact curves on M, locally equivalent to the solutions of a fourth-order 
ODE, are the geodesies of a sub-Riemannian metric only if a sequence of invariants vanish. 
The first of these, which was earlier identified by Fels, determines if the differential equation 
is variational. The next two determine if there is a well-defined metric on M and if the 
given paths are its geodesies. 



Introduction 

In this note, I will discuss the problem of recovering the geometric structure of a three- 
dimensional contact manifold with a sub-Riemannian metric from the geodesies for this 
metric. (Sub-Riemannian metrics are also known as Carnot- Carat heo dory metrics.) Since 
all the results herein will be local in nature, the manifold may be taken to be an open set 
U G R 3 with contact form dy — zdx, and we may assume that on contact planes the metric 
has the form 

g = Edx 2 + Fdxdz + Gdz 2 , 

where E, F, G are smooth functions on U such that g is positive definite. The geodesies, as 
constructed via the Griffiths formalism, form a collection of paths tangent to the contact 
structure, such that there is a 1-parameter family of distinct paths tangent to each contact 
direction at each point. Thus, part of the problem will be to determine which such collections 
of paths come from a sub-Riemannian metric. 

As explained below, the paths are locally equivalent to the integral curves of a scalar fourth- 
order ODE. The variational multiplier problem for fourth-order ODE — i.e., the problem 
of characterizing equations which are, up to multiple, the Euler-Lagrange equations for a 
second-order Lagrangian — was solved by M. Fels pi. Since sub-Riemannian geodesies arise 
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as solutions of a variational problem, the present work is an extension of that of Fels; to 
avoid confusion, the notation of || will be used whenever possible. 

I. Contact Path Geometries 

In this section I will review the construction of sub-Riemannian geodesies, and define a 
G-structure canonically associated to the geodesies as paths. 

Let M be an oriented three-manifold with contact distribution T> and sub-Riemannian 
metric g. It is standard that one can associate to g a SO (2)-structure N inside the oriented 
coframe bundle F{M), such that, for any coframing which is a local section of N, the forms 
(a? 1 , uj 2 , uj 3 ) of the coframing satisfy: 

(i) uj 3 annihilates the contact planes; 

(ii) {uj 1 ) 2 + {uj 2 ) 2 coincides with the metric on the contact planes; and, 

(iii) uj 1 A u 2 gives the induced orientation. 

Furthermore, one can choose N uniquely so that there is a connection form satisfying the 
structure equations^ 

du) 1 = (f) Au 2 + {aiuj 1 + a 2 uj 2 ) A uj 3 
duj 2 = -<p Auj 1 + {a 2 uj l - a x uj 2 ) A uj 3 
du) 3 = uj 1 A uj 2 
d(j) = Ku 1 A uj 2 mod uj 3 . 

The functions a 1; a 2 are components of the torsion of g and K is called its scalar curvature.0 
Every contact curve in M has a lift to iV on which the forms u 2 and us 3 vanish. (In this way, 
we'll identify N with the space of contact directions on M .) Applying the Griffiths formalism 
H to find the integral curves in N of the Pfaffian system {uj 2 ,uj 3 } which are extremal curves 
for arclength J uj 1 , we obtain the following characterization of sub-Riemannian geodesies^: 

Proposition 1.1. Let Z be the rank one affine subbundle ofT*N on which the canonical 
one-form is o = uj 1 — xuj 3 , x G M. (Forms on Z are pulled back via 7r : Z — > N .) Then 

1 This result appears in |^|, where it is attributed independently to Bryant-Hsu and to G. Wilkens. 

2 It's clear that taking the d 1 as an orthonormal coframe canonically associates to g with Riemannian 
metric g on M , which induces g on T>, and defines a canonical foliation perpendicular to T>. The torsion 
tensor is the Lie derivative of g along the leaves; if this vanishes, g descends to any (locally defined) quotient 
surface by foliation, and K is the Gauss curvature of the metric on that surface. 

3 See H for earlier derivations of the geodesies by other methods. 
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smooth geodesies are in 1-to-l correspondence, via the submersion Z — > M , with integral 
curves of the Pfaffian system T = {80, #1, 02, #3} on Z, where 

e = uj 3 

6, = ^ 

62 = <t> — xu 1 

63 = dx — a\U\ — O2002- 

Remark. In general, it is still an open question under what conditions all extremal curves for 
a given variational problem with differential constraints arise as projections of integral curves 
of the differential system formulated by Griffiths. For example, in sub-Riemannian geometry 
in dimension four, exceptional extremal curves exist which do not come from the Griffiths 
system. Essentially, this is because these abnormal minimizers || have few or no compactly 
supported variations that are tangent to the given distribution. However, by applying the 
regularity test given by Hsu f?| , one can show that for a sub-Riemannian metric on a contact 
manifold, all geodesies arise via the Griffiths formalism. (Intuitively, compact variations 
exist because contact curves can be locally expressed in terms of an arbitrary function and 
its derivatives.) 

Returning to the system T given above, let L be the line field on Z which is annihilated by 
T. Integral curves of this line field push down via 7r to give a 1-parameter family of curves 
through each point of N, and push down to M to give a 1-parameter family of geodesies 
tangent to each contact direction. We will now generalize this situation, throwing away the 
metric: 

Definition 1.2. Let M 3 be a contact manifold. Let p : P —>■ M be a fibration, with two- 
dimensional fibres, and L be a line field on P transverse to the two-dimensional fibres of p. 
Let I be the Pfaffian system on P which annihilates L and the fibres of p, and let J be the 
intersection of the retracting space 0] of I with the annihilator of L. Then (P,L,p) defines 
a contact path geometry on M if 

(i) for any vector v 6 L, p*(f ) is tangent to a contact plane; 

(ii) the first derived system X' is one- dimensional at each point of P; 
(Hi) J' = X at each point of P. 
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The last two conditions need explaining. Because of the transversality of L, X is two- 
dimensional. Because of condition (i), X' contains the pullback of any contact form on M. 
If X were integrable (i.e., X' = X, instead of being one-dimensional) then all paths through a 
given fibre p~ x {x) would project down to a single contact curve on M, so that there would 
be only one path through x G M. Condition (ii) implies that J is three-dimensional at each 
point. It is automatic that X C J' . If J were integrable, then integral surfaces of J would 
intersect p~ x {x) in a 1-parameter family of curves; since each such surface would project 
down to a single contact curve in M, this would imply that there was only a 1-parameter 
family of paths through x. 

Condition (ii) also implies that the three-dimensional distribution T> containing L and the 
kernel of p* is bracket-generating. That in turn guarantees, by Chow's theorem JT|, that two 
arbitrary points in M can be connected by a piecewise smooth sequence of paths. 

Proposition 1.3. Given a contact path geometry we can construct, in a neighbourhood of 
any point q G P , a coframe (a, O , 0i, 9 2 , 9 3 ) such that 

1. v G TP projects down to be a contact direction on M if and only if 9 (v) = 

2. X = {0 O) 0i} 

3. J = {0 O , 0i, 9 2 } 

4- L = {0 O , 0i, 2 , 03} 

Moreover, these forms satisfy 

(1) d9i = 9 i+1 A a mod O , . . . , 9 h < i < 2. 

These will be called 0- adapted coframes for the contact path geometry. 

Proof. Let p(q) = x G M. On a neighbourhood V of x, there exists a contact form 6q, and 

1-forms a, 9\ such that d9 = 9\ A a mod O . Pull these forms back to U — p~ l {V) C P; we 

will shrink U when necessary. Since a, 9\ both restrict to be zero along the fibres of p, they 

cannot be independent modulo X. Therefore we can arrange, by adding multiples of a, that 

0i G X. (Note that now 9\ is no longer the pullback of a form on V.). Since 9q, 9\ G L- 1 , 

then a £ L^. Since O G X', then d0i 7^ mod X. 

Since 0q, 0i, a span an integrable system, then there will be a smooth 1-form 2 on U that 

g?0i = 2 A a mod X. (Since X' is one-dimensional at each point, 02 is nonzero on U.) By 

adding multiples of a, we can arrange that 02 G L -1 , giving condition 3. Because J' 7^ ^7, 

4 



there must be a nonzero 1-form 63 on U such that d9 2 = 63 A a mod J . We can similarly 
arrange that #3 G L x . □ 



Remark. The above proposition could also be proved just using the assumption that L is a 
line field on P and X, ,7 satisfy conditions (ii,ifi) in Defn. |1.2| . The contact structure and 
the submersion to M can be recovered from X' and the retracting space C(X') respectively. 

Corollary 1.4. In some neighbourhood U of any given point q G P , there exist coordinates 
x, yo, yi, y 2 , 2/3 such that, for some function F on U, 

a = —dx 

9o = dy - yidx 

(2) 61= dyi - y 2 dx 

9 2 = dy 2 - y 3 dx 

O3 = dy 3 - F(x, yo, y±, y 2 , yz)dx 

is a 0-adapted coframe. Consequently, paths in P are locally equivalent to the solutions of 
the fourth- order ODE 

(3) y "" = F{x,y,y',y",y'"). 

Proof. The structure equations (0) enable us to apply the Goursat normal form theorem to 
system J . This gives J = {80, Q\, 9 2 }, in terms of the forms defined here. Since X = {6*0, #1}, 
then dx ^ L L , and so there exists some function F such that dy 3 — F(x, yo, yi, y 2 , y^)dx = 
along the paths in U. □ 

The set of 0-adapted coframes (cr,9o,0i,9 2 ,93) for given contact path geometry forms a 
principal bundle over P, with structure group Go C GL(5, R) consisting of matrices of the 
form 

(a * * \ 

6 

* a^b 

* * a~ 2 b 
^0 * * * a~ 3 b J 

This is precisely the G-structure that Fels associates to a fourth-order ODE up to contact 

transformation (cf. Lemma 3.1). Since the path geometry can be recovered uniquely 

from the G-structure, we will treat the two notions as synonymous. 
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2. Variational and Sub-Riemannian Path Geometries 



The goal of Cartan's method of equivalence is, for a given G-structure, to find a sub- 
bundle, with reduced structure group, on which there exists a unique connection. Like the 
Levi-Civita connection in Riemannian geometry, this is typically obtained by fixing the value 
of all or part of the torsion of the connection. Then, invariants may be extracted from the 
remaining torsion or the curvature of the connection. 

We begin with Fels' result for Go-structures of coframes satisfying (|1|). This gives a 
reduction of structure to the subgroup G\ C Go consisting of matrices of the form 



(a 
b 








a~ x b 
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In terms of path geometry, the result is: 
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Theorem 1 (Fels ||). Let Bq \ P define a contact path geometry. Then there is a sub- 
bundle E>i with three-dimensional structure group G\, on which there exists a unique equi- 
variant connection satisfying the following structure equations: 

e 1 A{T 4 e 2 + T s e 3j 



(4a) da = a A a + O A (TA + T 2 2 + T 3 9 3 ) 

(4b) d0 = (3 A 6 + a A 6 1 

(4c) dB\ = {(3 - a) A 0i + 7 A O + a A 6 2 

(4d) d0 2 = ((3 - 2a) A 6 2 + §7 A B x + o A 6 3 

(4e) d6 3 = (13 - 3a) A 6 3 + 7 A 6 2 + a A (J O O + ^i) + T 6 6 A X + T 7 o A 2 + T 8 0! A 6 2 . 

[The one-forms a, /3, 7 are connection forms, and Iq, Ji, Ti, . . . T 8 are components of the tor- 
sion of the connection.] 

Moreover, assuming P is locally defined by a fourth-order ODE (|3|) ; solutions of that ODE 
are critical curves for a second-order Lagrangian if and only if the relative invariants Ii and 
T 5 both vanish identically on B\. In that case, T 8 also vanishes. 



The essence of Fels' proof of the second statement is exhibiting a two-form on B\, 



u = m(e A9 3 - 0i A0 2 ), 
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where m is a non-zero function, such that uj is closed and Gi-invariant. (In fact, dlogm = 
3a — 2/3, and, as Fels notes, the structure equations imply that that one-form is closed in 
the variational case.) It then follows that u is the exterior derivative of the Poincare-Cartan 
form associated to a Lagrangian on the space of 2-jets. 

We will speak of a path geometry for which ii, T 5 , T 8 vanish identically as being variational. 

Example 1. Consider the second-order Lagrangian J e~ 3y "dx, for which the Euler-Lagrangian 
equations are, up to multiple, 

y"" - 3(y"') 2 = 0. 

The coframe (Q) gives a section of the bundle B defining the corresponding G -structure 
on J 3 (M, M). This coframe may be modified to give the following section of the reduced 
structure B\. 

do = dy - yidx 
#i = dyi — y 2 dx 
(5) 2 = dy 2 - y 3 dx - y 3 9 1 + ^y^o 

03 = dy 3 - 3y 3 dy 2 - ±y$6i + §2/f0 o 
o = dx + 0i - §2/ 3 6>o 

Of course, the torsion satisfies I\ = T 5 = T 8 = 0, but one also may computeQthat T 2 = ^-y 3 , 
T 3 = | and T 4 = — 1 along this section of B\. 



Example 2. (sub-Riemannian geometry) Let Z be the five-manifold of Proposition It 
is easy to verify that the 1-forms given there, when rounded out by a = uj 1 — xco 3 , form a 
0-adapted coframe for the corresponding contact path geometry. We may adapt the coframe 



4 In order to evaluate the torsion components along a given section of B\ , one must determine the values 
of the connection forms in terms of the given coframe. To do this, begin with the d0Q equation (4b), which 



determines (3 modulo 8q. One may set (3 = /3q + b0Q, where /3o is any form satisfying (ph|) and b is not yet 
determined. Then ( fig 1 ) determines a and 7 modulo 0q,9i. In fact, one may set 

a = da + aOa + zB\ 
7 = 70 - a9i + c9 . 



Now ([la]) determines z while (^d|) , (ffe|) determine a, b and c. 
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to obtain a section of the reduced bundle B\ \ Z: 

9 = u 3 

e l = u 2 i 

A = — (a 2 + 3x 2 - 3K) 

(6) 9 2 = <f ) -xu 1 + Auj 3 with 10 

3 = dx- ai uj l - (a 2 + A)uj 2 + Bcu 3 B = — (s 2 - 3fci - 6aiX - 216 x ) , 

1 3 o 

a = to xoo 

5 

where K is the scalar curvature, and the b iy Sj and fcj are defined on iV by 
dai = 2a 2 (f) + (si + 6 2 )^ 1 + (s2 + &i)^ 2 
(ia 2 = — 2ai0 + (s 2 — + (b 2 — Si)lu 2 } mod u 3 
dK = kiuo 1 + k 2 uj 2 

Again, one may compute that I\ = T5 = Tg = 0, confirming that the path geometry is 
variational, while T 2 = 0, T 3 = |, and T 4 = —1 for this coframe. 

The fact that we obtained the same values for T 3 and T 4 as those from a general second- 
order Lagrangian hints at further relations among the torsion components. One uncovers 
one of these while deriving the refined structure equations: 

Proposition 2.1. Let B\ be the canonical G\-structure for a variational path geometry. 
Then there exist functions Ui,U 2 on B\ such that the connection forms satisfy 

da = \d3 

3 r = T\Q\ + T 2 9 2 + T 3 9 3 

(7) d3 = a A 7 - r A 9 1 - 3v A 9 where 

v = U l 9 1 + U 2 9 2 - T 2 9 3 + T 7 a 

(1) e 7 A a - r A0 2 — v A9 1 mod 9 
The torsion components satisfy T 3 = — |T 4 and 

(8) dT x = 71 (2a - 2(3) - §T 27 - 2U x o mod 9 , 6 U 9 2 , 9 3 

(9) dT 2 = T 2 (3a - 23) - §T 47 - (7\ + 2U 2 )a mod O , 0i, fc, 9 3 

(10) cff 4 ee T 4 (4a - 23) - §T 2 cr mod 9 , 6 U 9 2 . 

The above equations indicate that T 4 is a relative invariant on i.e., it varies along the 
fibres only by scaling. Moreover, they indicate that the quadratic form g = a 2 — T A 9\ is 



well-defined, up to multiple and modulo 0o, on N. For, suppose v is a vector field on Bi 
which is annihilated by a, 9q, 0i, 9 2 - Then computing the Lie derivative of g gives 

C v {g) = 2a o ( v da) - 2T 4 0i o (v -. d^) - (t» -. dT 4 )0 2 
= 2(v -«a) [a 2 - T 4 2 ] mod O . 

This quadratic form will be our candidate for a sub-Riemannian metric. Matters being so, 
we will say that a variational geometry is nondegenerate if T 4 7^ everywhere,^ and definite 
if T 4 is negative everywhere. Assuming the latter is the case, then we may normalize T 2 and 
T 4 to have the same values as in Example [|. 

Proposition 2.2. Let B\ \ P be a definite variational path geometry. Then there is a 
sub-bundle B2 C B\ on which 

T 2 = and T 4 = — 1 . 
On B 2 there exist smooth functions W , Wi, W 2 , G ,Gi,G 2 ,G 3 , and H such that 

(11) (3 = 2a + W 9 + Wi0i + W 2 9 2 

(12) j = Ha- 3(G 9 + + G 2 9 2 + G 3 9 3 ). 

Proof. Structure equations (PD,(|I0D show that we may first pass to the sub-bundle where 
T 4 = — 1 and then move along the fibres in a direction dual to 7 to pass to the sub-bundle 
where T 2 = 0. Once there, these equations show that (3 — 2a and 7 restrict to have the above 
form. Of course, @ shows that |if = T\ + 2U 2 . □ 

On B 2 , the structure equations (f|) take the form 

da = a A a + 9 A (TA + §0 3 ) -9 t A9 2 
d9 = 2a A O + cr A 6>i 

(13) d0i = (/3-a) A0i+7A0 O + ^A0 2 
d9 2 = (/3 - 2a) A 2 + §7 A 9 1 + a A #3 

d0 3 = (/? - 3a) A 3 + 7 A 2 + /off A O + T 6 o A 9 X + T 7 o A 2 . 
with given by ([TT|) . 

5 Suppose a variational path structure has T3 and T 4 identically zero; the refined structure equations show 
that T2 = also. Recall that the system which restricts to be zero along the fibres of p : P — > M 3 is spanned 
by <7, 9o, 9i- Since dcr = Ti^o A 9\ mod cr, vectors that are in the kernel of a push down to give a well-defined 
plane field on M. These planes intersect the contact planes in a distinguished family of contact directions, 
which are null lines with respect to g. 
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It's clear that the fibres of B 2 are one-dimensional, and the structure group of B 2 is simply 
K*. A element A 7^ of this group acts on sections of B 2 by 

9\ ■ (a, , 0i, ^2, 0a) = (Ac, A 2 o , A0 X , 6 2 , X~%). 

Structure equations (0) show that g\a = a, g* x [3 = j3, and 5^(7) = A _1 7. Then the action on 
the new torsion is clearly 

9x ■ (W , W x , W 2 ) = (\~ 2 W , X-'Wu W 2 ), 
gx • (H, Gq, G±, G 2 , G3) = (A 2 H, A 3 G*o,A 2 Gi,A 1 G 2 ,Gs). 

In particular, W 2 , G 3 , and the ratios G\ : Wq and G 2 : H 7 ! are invariant under the scaling 
action. 

We should expect this scaling to be present, since two sub-Riemannian metrics which differ 
by a constant factor have the same geodesies and hence define the same path geometry. For 
purposes of constructing a specific metric, we will need to choose a section of B 2 . Since 
3a — 2/3 is closed, integrals of this one-form comprise a canonical codimension-one foliation 
of B 2 which is transverse to fibres and invariant under the scaling action. 

Definition 2.3. A section of B 2 along which 

(14) 3a - 2(3 = 

will be called a canonical section of B 2 , or a canonical coframe on P. It follows from fllTD 
that 

(15) a = -2(I^o0o + Wi0i + ^ 2 2 ) 
along a canonical section. 

One can check that the coframing constructed in Example ^| is a canonical coframe. Since 
such coframings are unique up to scale, it follows that if a path geometry comes from a sub- 
Riemannian metric, then in terms of a canonical coframe that metric must be g = a 2 + (0i) 2 . 

Proposition 2.4. Let P be a definite variational path geometry on contact manifold M 3 
and (<t,0o,0i,02,03) a canonical coframe on P. Then g = a 2 + (0i) 2 gives a well-defined 
metric on the contact planes of M if and only ifW 2 is identically zero on P. 
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Proof. Let v be any vector field on P tangent to the fibres of the projection p : P —>■ M. 
Since v is annihilated by Oq, 0\ and a, 

C v {g) = {v^a){af + {v^{(3-a)){d l f mod 6 
= -W 2 (v^9 2 ) (2(a) 2 + {O,) 2 ) . 

□ 

Although the coframe (|5|) is not a section of B 2 , it can be adjusted so that T 2 = 0, 
whereupon we see that W 2 is nonzero for Example |I|. 

For the rest of this section we will assume that W 2 is identically zero. It remains to be 
seen if the given paths on M 3 — which are projections of the integral curves of the line field 
L — are the geodesies of the sub-Riemannian metric we have constructed. To investigate this 
further, we will need the torsion identities 

G 3 = 0, G 2 = W h 

which result from computing d(d9i) = using the structure equations fll3|) and the equations 



with W 2 = 0. 

Remark. One might wonder if other identities hold among the remaining torsion coefficients 
Go, Gi, H, I , Ti, Tq, T 7 , Wo, W\ as a result of our assumption that W 2 = 0. However, no 
further identities arise, and this is proved by showing that the exterior differential system 
defining a G-structure satisfying the structure equations on B 2 with W 2 = is involutive. 

Theorem 2. Let P carry a variational and definite path geometry with invariant W 2 iden- 
tically zero, and let {o~,6o,0x,9 2 ,6z) be a fixed canonical coframe on P. Then the paths in 
P project to be geodesies in M for the sub-Riemannian metric of Prop. \2.4\ if and only if 
G\ = 2Wo identically on P. 

Proof. Let N be the quotient of P by the foliation by integral curves of the system = 
{a, 9o, 0i, 6 2 }. Each contact curve in M has a unique lift to iV as an integral curve of 
X = {#o, Clearly, arclength is measured along these lifts by the integral of a modulo X. 
However, the form a on P does not descend to be well-defined on N, as shown by 

da = a A a + 6 A (TA + T 3 9 3 ) + T A Q X A Q 2 

= §0 O A # 3 mod A 2 X (1) . 
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However, computing d 2 9 = shows that 

(16) dW x = (Gi - W )a + mod X, 
and this, together with d9 = mod A 2 I^\ shows that the 1-form 

a = a + 3Wi0 o 

is well-defined on N. 

Now arclength with respect to the metric may be measured on the integral curves of X by 
the Lagrangian J a. We will apply the Griffiths formalism || to investigate which of these 
are geodesies for g. Then, we will try to find conditions under which these curves coincide 
with the projections of the paths in P under 7r : P — > N. 

Let £ = a + x9q + yQ\ on Y = N x R 2 . Then one finds that the two-form d^ is of full rank 
on Y, except where y = 0. Accordingly, let £ = a + a;^ on Z = JV x i. Now one computes 
that 

(17) d£ = (dx + (3d - Wo)<r) A 6 + (0 2 + (x + SWJa) A ^ mod 6 A ^. 

Let /C be the rank four Pfaffian system on Z spanned by the four one-forms on the right in 
(0): 

K. = {#o, u 9 2 + (x + 5Wi)<7, rfx + (3Gi - l^o)^}- 
According to the Griffiths formalism, integral curves of /C project to be extremal curves 
for J a on N. These coincide with the projections of the paths in P if and only if, in a 
neighbourhood U of each point of P, there is a local diffeomorphism ip : U — > Z such that 
<^*/C coincides with L -1 = {# , ^2, ^3}, the Pfaffian system on P which defines the paths. 
(The diffeomorphism would follow from the identification of paths with geodesies on N.) 
The form (p*{6 2 + (x + 5Wi)cr) belongs in L x if and only if ip*x = —5W\. Then, by (|16D , 

<p*(dx + (3d - W )<t) = (4W - 2Ga)o- mod L x , 

showing that v?*/C = L 1 - if and only if G x = 2W . □ 

3. Discussion 

The results of the previous section may be surprising. For, one could reason that, once a 
path geometry is known to be variational, it must arise from a second-order Lagrangian, of 
the form 

J L(x,y,y',y")dx, 
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satisfying the nondegeneracy condition d 2 L/d(y") 2 ^ 0. Then L is of the form 

L(x,y,y',y") = \J E + Fy" + G(y") 2 , 

for some functions E, F, G of x, y, y', if and only if L satisfies a certain third-order ODE as 
a function of y" . In other words, it seems like only one extra condition must be satisfied in 
order for the Lagrangian to be length with respect to a sub-Riemannian metric. Instead, 
we find that two scalar conditions (in addition to the Fels variational condition) must hold 
in order for the metric to be well-defined and in order for its extremals to coincide with 
the given paths. (The reader should note that the above remark about involutivity implies 
that the condition G\ = 2Wq is independent from W2 = 0.) It would be interesting to 
find examples of variational path geometries (equivalently, variational fourth-order ODE) 
for which W% = but the extremals of the associated metric do not coincide with the given 
paths. Such examples must exist, again, because of the involutivity of W2 = 0. 
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